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ABSTRACT 

We present new L' (3.75 /im) photometry of six L and T dwarfs, and M' (4.70 /im) photometry 
of ten L and T dwarfs, observed at Gemini (North) Observatory, and new 3.55, 4.49, 5.73 and 
7.87 /im photometry of nine L and T dwarfs, obtained with the Spitzer Space Telescope. The 
sample includes unusually blue and red dwarfs selected from our near-infrared studies. The data 
are combined with previously published L' , M' and Spitzer photometry of L and T dwarfs, and 
trends of colors with spectral type and other colors are examined. Recent model atmospheres 
by Marley and Saumon are used to generate synthetic colors for ranges of effective temperature, 
gravity, grain sedimentation efficiency, metallicity and vertical mixing efficiency. We explore 
how these parameters affect the mid-infrared colors of L and T dwarfs and find that the data 
are modelled satisfactorily only if substantial vertical mixing occurs in both L- and T-dwarf 
atmospheres. The location and range of the L and T dwarf sequences in IRAC color-color and 
color-magnitude diagrams is also only reproduced if this mixing occurs, with a range of efficiency 
described by K zz ~ 10 2 -10 6 cm 2 s _1 . The colors of the unusually red dwarfs are best reproduced 
by non-equilibrium models with low sedimentation efficiency, i.e. thick cloud decks, and the colors 
of the unusually blue dwarfs by non-equilibrium models with high sedimentation efficiency, i.e. 
thin cloud decks. The K—L' and Spitzer [3.55]-[4.49] colors can be used as indicators of effective 
temperature for L and T dwarfs, but care must be taken to include gravity and metallicity effects 
for late-T dwarfs and vertical mixing for both lgte-L and T dwarfs. 
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1. Introduction 

The Sloan Digital Sky Survey (SDSS; York et 
al. 2000) and the Two Micron All Sky Survey 
(2MASS; Bcichman et al. 1998, Skrutskie et al. 
2006) have revealed large numbers of ultracool 
low-mass field dwarfs, known as L and T dwarfs. 
The effective temperatures (T e ff) of L dwarfs are 
~ 1450-2200 K, and those of currently known 
T dwarfs are - 700-1450 K (Golimowski et al. 
2004a, hereafter G04; Vrba et al. 2004). 

As T c g drops below that of the late-M dwarfs, 
the spectral energy distributions (SEDs) of ultra- 
cool dwarfs are dramatically changed by chemical 
and cloud processes in their atmospheres. First, 
iron and silicate grains condense high in the at- 
mosphere and form clouds that veil gaseous ab- 
sorption bands and redden the near-infrared colors 
of L dwarfs. As T e s falls, these clouds form pro- 
gressively deeper in the atmosphere and become 
more optically thick. By ~ 1500-1700 K, the ef- 
fect of the clouds on the emergent flux is great- 
est (Ackerman & Marley 2001). At lower T eff , the 
clouds lie near or below the base of the wavelength- 
dependent photosphere and only marginally affect 
the flux distributions of T dwarfs. The rapidity, 
with respect to T ff, with which the clouds di- 
minish is not well understood and several mech- 
anisms have been proposed (e.g. Burgasser et 
al. 2002b, Burrows, Sudarsky & Hubeny 2006, 
Knapp et al. 2004, hereafter K04). G04 showed 
that T c ff is nearly constant (<~ 1300-1500 K) for 
spectral types L7 to T4. Thus, for these spectral 
types, the primary cause of the observed spec- 
tral change is apparently rapid alteration of the 
vertical properties and distribution of the cloud, 
over a small range in T c ff. The changing cloud 
properties alter the optical depths and brightness 
temperatures from which the flux emerges. As 
the clouds depart, the upper atmosphere cools 
and CH4 supplants CO as the dominant carbon- 
bearing molecule. CH4 appears in the 3-4 /zm 
spectra of mid-L dwarfs (Noll et al. 2000), the 2.2- 
2.5 and 8 /im spectra of late-L dwarfs (Geballe et 
al. 2002; Roellig et al. 2004; Cushing et al. 2006), 
and the 1.6 fim spectra of TO dwarfs (Geballe et 
al. 2002). Together, increasing CH4 absorption 
and sinking cloud decks cause progressively bluer 
near-infrared colors of T dwarfs. 

Over 500 ultracool field dwarfs have been con- 



firmed and studied in the near-infrared (e.g. Bur- 
gasser et al. 2006b, Chiu et al. 2006 (hereafter 
C06), Kendall et al. 2004, Tinney et al. 2005). 
Our group has so far obtained accurate J, H, and 
K photometry on the Mauna Kea Observatories 
(MKO) system (Simons k Tokunaga 2002, Toku- 
naga et al. 2002) and R w 500, 0.9-2.5 fim spec- 
troscopy for 102 L dwarfs and 65 T dwarfs. This 
sample is summarized by K04 and C06. Leggett 
et al. (2002) and G04 extended our study to 
longer wavelengths, presenting L' (3.75 /xm) and 
M' (4.70 ^m) photometry for a subset of the JHK 
sample. In this paper, we extend our study further 
by presenting broadband photometry for wave- 
lengths up to 7.9 ^m. 

The 3-15 ^m mid- infrared region is interest- 
ing for ultracool-dwarf studies for many reasons. 
This region contains molecular absorption bands 
of CH 4 , CO, H 2 and NH 3 (Cushing et al. 2006) 
that are prominent in the spectra of ultracool 
dwarfs at various temperatures. Figure 1 shows 
observed and synthetic 3-10 /im spectra for late- 
ly and T dwarfs and identifies the principal ab- 
sorbing species. The vibrotational line lists for 
CH4, NH3 and H2O are more complete in the 
mid-infrared than the near-infrared, so the uncer- 
tainties in the synthetic mid-infrared spectra are 
smaller than those in the near-infrared spectra. 
Also, the mid-infrared flux emerges from photo- 
spheric regions that have a much smaller range of 
brightness temperatures than those giving rise to 
the near-infrared flux (Ackerman & Marley 2001). 
Hence the mid-infrared spectrum should be easier 
to interpret than the near-infrared spectrum. Fur- 
thermore, mid-infrared photometry is more sen- 
sitive than near-infrared photometry to vertical 
mixing of trace species, which can result in pho- 
tospheric chemical abundances very different from 
equilibrium (Saumon et al. 2003; G04). 

Obtaining mid-infrared photometry and spec- 
troscopy with ground-based telescopes is difficult, 
due to the high and rapidly varying background at 
these wavelengths, and to the strong telluric ab- 
sorption bands of CH4, CO, CO2, H2O, N2O and 
O3. The Spitzer Space Telescope (Werner et al. 
2003) has made possible observations at these and 
longer wavelengths that were previously impossi- 
ble from the ground. We are using Spitzer's In- 
frared Array Camera (IRAC; Fazio et al. 2004) and 
Infrared Spectrograph (IRS; Houck et al. 2004) 
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to obtain 3.6-7.9 /jm photometry and 6-15 /xm 
low-resolution spectra of late-L and T dwarfs. In 
this paper, we present the IRAC results along with 
ground-based V and M' photometry obtained us- 
ing the Near InfraRed Imager and spectrograph 
(NIRI; Hodapp et al. 2003) at the Gemini (North) 
Observatory. The IRS spectra will be presented in 
a future paper. 

2. Observations 

2.1. Gemini NIRI L' and M' Photometry 

We obtained MKO L' and/or M' images of 12 
late-L and T dwarfs, listed in Table 1, on 10 dry 
and photometric nights between UT 2004 April 05 
and UT 2005 January 28. The data were obtained 
in queue-observing mode through programs GN- 
2004A-Q-16 and GN-2004B-Q-43, using Gemini 
Observatory's NIRI. The half-power bandpasses 
for the MKO V and M' filters (Simons & Toku- 
naga 2002, Tokunaga et al. 2002) are shown in 
Figure 1 and listed in Table 2. Because the sky is 
bright at these wavelengths, we used NIRFs small- 
est pixel scale (0'.'022 pix -1 ) to avoid saturating 
the detector. The full array readout was used, pro- 
viding a field of view of 22'.'5 x 22'.'5. Given the 
small field of view, we acquired each target using 
proper-motion-corrected coordinates (if available) 
and a short J-band verification image. During 
the L' and M 1 exposures, the target was sequen- 
tially offset from the center of the detector toward 
each of the four corners, generating a four-position 
dither pattern with 6" offsets. Offset guide stars 
were used. Each L' image comprised 24 coadded 
exposures of 0.8 s, and each M' image comprised 
34 coadded exposures of 0.5 s. The total observa- 
tion time per target was typically 1.5 hr, most of 
which was spent on the M' exposures. Each set 
of science images was accompanied by contempo- 
raneous dark frames and images of bright photo- 
metric standard stars (Leggett et al. 2003) taken 
with the same filters. Flat fields were generated 
by median-stacking the dithered science images. 

The data were reduced using a version of the 
astronomical imaging pipeline ORAC-DR (Ca- 
vanagh et al. 2003) modified for NIRI. 1 The 

1 The ORAC-DR Imaging Data Re- 

duction User Guide is available at 
littp: / /www. starlink.ac.uk/star/docs/sun232.htx/sun232. html. 
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Fig. 1. — Observed and synthetic spectra of late- 
L and T dwarfs (normalized to peak flux) in the 
wavelength interval where we report new photom- 
etry. The black curves are NIRI and IRS spec- 
tra of the L9 dwarf 2MASS J09083803+5032088 
(Stephens private communication). The red curve 
is a synthetic spectrum for T c s = 1400 K, / sc d = 2, 
and log g — 5, representing late-L types, and the 
blue curve is a synthetic spectrum for T c g =1000, 
/scd = 4, and log 5 = 5, representing mid-T 
types (Marley & Saumon, private communica- 
tion). Strong molecular absorption bands are indi- 
cated, as are the bandpasses of the MKO L' , MKO 
M' and IRAC channels 1-4 filters: [3.55], [4.49], 
[5.73], [7.87]. 
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pipeline uses Starlink routines to perform bad- 
pixel masking, dark subtraction, differencing of 
adjacent pairs, flat-field creation and division, 
feature detection and matching between object 
frames, and resampling. The final combined im- 
age, termed a mosaic, is generated with intensity 
offsets applied to give the most consistent results 
in the overlapping regions of the dithered images. 
Pairwise differencing allows more accurate sky 
subtraction and flat-fielding, although it yields 
positive and negative images of each source in the 
mosaics. The mosaics were binned 3 x 3 to en- 
hance the visibility of the L and T dwarfs. The 
final pixel scale was therefore 0'.'066 pix^ 1 . 

We performed aperture photometry on both the 
positive and negative images of the L and T dwarfs 
using apertures of diameter l"-2". We defined 
the uncertainty to be the difference between the 
measurements of the positive and negative images. 
This difference is similar to, or slightly larger than, 
the formal error calculated from the noise. Table 1 
lists the resulting V and M' magnitudes, as well as 
any previously published measurements. Hence- 
forth, we use the weighted averages of the new 
and previously published measurements of all the 
dwarfs except 2MASS J15031961+2525196, whose 
previous uncertain M' measurement is discarded 
in favor of the NIFJ value. 

2.2. Spitzer Four-Channel Photometry 

We obtained IRAC four-channel photometry 
(3.55, 4.49, 5.73, and 7.87 fj,m) of the nine L and T 
dwarfs listed in Table 3 between UT 2005 July 23 
and UT 2006 July 06. The data were obtained as 
part of our Cycle 2 Spitzer General Observer pro- 
gram # 20514. The half-power bandpasses for the 
four channels are listed in Table 2 and indicated 
in Figure 1. All four channels have 256 x 256-pixel 
detectors with a pixel size of l'.'2xl'.'2, yielding a 
5.'2x5.'2 field of view. Two adjacent fields are im- 
aged in pairs (channels 1 and 3; channels 2 and 
4) using dichroic beam splitters. The telescope is 
then nodded to image a target in all four channels. 
2 We used exposure times of 30 s and a 5-position 
medium-sized (52 pixels) dither pattern repeated 
one to four times. The total observing time per 
target ranged from 11 to 26 minutes. The full ar- 

2 See Fazio et al. (2004) and the IRAC Users Manual at 
http://ssc.spitzer.caltech.edu/irac/descrip.html 



ray was read out. 

The data were reduced using the post-basic- 
calibration data mosaics generated by version 14 
of the IRAC pipeline. 3 The mosaics were flat- 
fielded and flux-calibrated using super-flats and 
global primary and secondary standards observed 
by Spitzer. We performed aperture photometry 
on the L and T dwarfs using apertures of radii 
3-7 pixels. The dwarfs were usually isolated and 
well detected, so apertures of 7 pixels could be 
used. We applied channel-dependent aperture cor- 
rections as described in Chapter 5 of the IRAC 
Data Handbook. 3 For channels 1 to 3, the correc- 
tion for the 7-pixel aperture is 3%, for the 5-pixel 
aperture it is 5-6% and for the 3-pixel aperture 
it is 11-12%. For channel 4, these corrections 
are 4%, 7%, and 22%. We derived photometric 
errors from the uncertainty images that are pro- 
vided with the post-basic-calibration data. The 
photometry was converted from milliJanskys to 
magnitudes on the Vega system using the zero- 
magnitude fluxes given in the IRAC Data Hand- 
book (280.9, 179.7, 115.0, and 64.1 Jy for channels 
1 to 4, respectively). The fluxes, magnitudes, and 
errors for the observed dwarfs are given in Table 3. 
Note that in addition to the errors in Table 3, 
there are absolute calibration uncertainties of 2- 
3%. There are also similarly-sized systematic un- 
certainties introduced by pipeline dependencies, as 
we found by comparing the photometry produced 
by pipeline versions 12, 13 and 14 for seven of our 
targets. Photometry produced by versions 12 and 
14 agreed well, to typically 4%, however version 13 
produced fainter magnitudes than the other two 
pipelines, by ~7%. The Spitzer helpdesk team 
recommend version 14 of the pipeline, and sug- 
gest that version 13 used an over-aggressive pixel 
rejection algorithm. We adopt the total photo- 
metric uncertainty to be the sum in quadrature of 
the values given in Table 3 plus 3%. 

IRAC flux densities follow the convention es- 
tablished by the Infrared Astronomical Satellite 
(IRAS) and other missions. The IRAC flux den- 
sity at a nominal wavelength assumes that the 
source has a spectrum described by F v oc v . 
However, this assumption is not valid for the SEDs 
of latc-L and T dwarfs (Figure 1); A x F\ is not 

3 Further information can be found at 
http:/ /ssc. spitzer. caltcch.edu/irac/dh/ 
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constant in these cases. Hence the fluxes given 
in Table 3 are not representative of the nominal 
wavelengths given in Table 2. For example, the 
3.55 /xm and 7.87 ^m IRAC fluxes of the L9 dwarf 
shown in Figure 1 differ by ~ 30% and ~ 10%, 
respectively, from the actual fluxes at those wave- 
lengths measured by Spitzer's IRS. If the source 
spectrum is known, the quoted IRAC flux for a 
nominal wavelength Aq can be expressed as 



flKAV (\ 



»> = / 



(y Q /v)f{y)S{y) dv / / {v Q /v) 2 S{v) dv 



where f(v) is the source spectrum and S{v) is the 
system response function (Cushing et al. 2006). 

To put our data into context, we augment our 
sample with the IRAC photometry of M, L, and 
T dwarfs from Patten et al. (2006, hereafter P06), 
and include the photometry of two T dwarfs dis- 
covered in IRAC images by Luhman et al. (2006). 

3. Colors and Spectral Types 

Figure 2 shows the differences between IRAC 
channel 1 (3.55 /zm) and MKO L' (3.75 ^m) mag- 
nitudes and between IRAC channel 2 (4.49 /im) 
and MKO M' (4.70 /im) magnitudes, as functions 
of spectral type. The data are from this work, 
Leggett et al. (2002), Reid & Cruz (2002), G04, 
and P06. IRAC channels 1 and 2 span wavelengths 
similar to those of L' and M', respectively (Fig- 
ure 1), and correlations can be used to check the 
consistency of Spitzer and ground-based data. 

Figure 2 shows that the IRAC channel 1 mag- 
nitudes (denoted [3.55]) become fainter than L' 
magnitudes for later spectral types. This behav- 
ior is expected because channel 1 samples more 
of the strong center of the 3.3 /xm CH4 absorption 
band, and less of its long- wavelength wing and flux 
peak near 4 /im, compared to the L' bandpass (see 
Figure 1). The T7.5 dwarf 2MASS J12171110- 
0311131 is clearly anomalous in this plot; we dis- 
cuss it further in §4. The difference between 
IRAC channel 2 ([4.49]) and M' is fairly constant 
([4.49]-M' « -0.3) for all types, with the ex- 
ception of the extremely red late-L dwarf 2MASS 
J22443167+2043433 (§4). Both channel 2 and M 1 
sample H2O and CO bands, even among T dwarfs 
(§4). Although channel 2 is much wider than M', 
the filters are similarly sensitive to changes with 
spectral type. 
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Fig. 2.— Comparison of MKO V (3.75 /zm) 
and M' (4.70 fim) magnitudes with the IRAC 
channel 1 ([3.55]) and channel 2 ([4.49]) magni- 
tudes. Triangles represent L dwarfs and squares 
represent T dwarfs. Filled symbols are data 
presented in this work; open symbols are data 
from Leggett et al. (2002), Reid & Cruz (2002), 
G04, and P06. The following dwarfs are identi- 
fied with abbreviated names: 2MASS J12171110- 
0311131, 2MASS J1534498-295227AB and 2MASS 
J2244316+204343. 
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Fig. 3. — Colors as a function of spectral type. The symbols are described in Figure 2. Here, and in the fol- 
lowing figures, open symbols are data from Leggett et al. (2002), Reid & Cruz (2002), G04, P06 and Luhman 
ct al. (2006). The following unusually blue or red dwarfs are identified and discussed in the text: 2MASS 
J0937347+293142, 2MASS J12171110-0311131, 2MASS J2244316+204343, SDSS J080531. 80+481233, SDSS 
J133148. 90-011651.4, SDSS J151643.01+305344.4. 



6 



Figure 3 shows various IRAC and MKO KL'M' 
colors plotted as functions of spectral type. The 
data are taken from this work, Leggett et al. 
(2002), Reid & Cruz (2002), K04, G04, C06, Luh- 
man et al. (2006) and P06. All of the colors show 
scatter at the ^0.5 mag level from mid-L to early- 
T types, but general trends and some peculiar- 
ities can be seen. 2MASS J22443167+2043433, 
the reddest L dwarf currently known (Dahn et 
al. 2002; K04), is clearly anomalous in many 
colors. SDSS J080531.80+481233.0 and SDSS 
J133148.90-011651.4, both unusually blue in the 
near-infrared (K04), are also anomalously blue 
in the IRAC [4.49]-[5.73] and [3.55]-[4.49] colors. 
SDSS J151643.01+305344.4, a red early-T dwarf 
(C06), is also red in #-[3.55] and [5.73] -[7.87]. 
Thus, the unusual traits seen in these dwarfs at 
shorter wavelengths often persist at longer wave- 
lengths. 2MASS J09373487+2931409, which is ex- 
tremely red in K-L' and #-[3.55], has suppressed 
if-band flux most likely due to very strong H 2 ab- 
sorption (Burgasser et al. 2002a). 

4. Interpretation of Colors and Luminosi- 
ties 

In this section, we compare the observed IRAC 
and MKO colors with the predicted colors from 
model atmospheres, and investigate the physical 
parameters that can produce the observed trends 
and variations seen in Figure 3. 

4.1. Description of Model Atmospheres 

State-of-the-art models of ultracool-dwarf at- 
mospheres have been developed by members of 
our team (Ackcrman & Mar ley 2001; Marley et 
al. 2002; Saumon et al. 2003; and will be fur- 
thered described in an upcoming paper). These 
models yield temperature-pressure-composition 
structures under conditions of radiative-convective 
equilibrium using the thermal radiative transfer 
source function technique of Toon et al. (1989). 
The gas opacity database includes the molecular 
lines of H 2 0, CH 4 , CO, NH 3 , H 2 S, PH 3 , TiO, VO, 
CrH, FeH, C0 2 , HCN, C 2 H 2 , C 2 H 4 , and C 2 H 6 , 
the atomic lines of the alkali metals (Li, Na, K, 
Rb and Cs), and continuum opacity from H 2 col- 
lisionally induced absorption (CIA), H 2 , H and 
He Rayleigh scattering, bound-free opacity from 
H and , and free- free opacity from He, H^ , and 



H2 • The models also allow opacity from arbitrary 
Mie-scattering particles. 

The radiative transfer code is closely coupled to 
cloud models that use a sedimentation efficiency 
parameter, / sc( j, to describe the balance between 
the downward sedimentation of condensates and 
the upward convection that replenishes the grains. 
Large values of / sc d correspond to rapid growth 
of large grains, which then quickly fall out of the 
atmosphere and yield physically and optically thin 
clouds. When / sc( j is small, the grains grow more 
slowly and the clouds are thicker. The models 
have been successfully applied to the atmospheres 
of brown dwarfs and Jupiter (Ackerman & Marley 
2001; K04; G04). 

The models have been further advanced to in- 
clude vertical mixing, which draws up chemical 
species from the well-mixed, deep, hot layers into 
the cool radiative photosphere. Although the pho- 
tosphere would be expected to be stable and in 
chemical equilibrium, even a small degree of mix- 
ing (such as is commonly observed in the strato- 
spheres of the solar system's giant planets) en- 
hances the abundances of stable species like CO 
and N 2 relative to their chemical equilibrium val- 
ues, while those of H 2 0, CH4 and NH3 are reduced 
(Fegley & Lodders 1996, Lodders & Fegley 2002, 
Saumon et al. 2003, 2006). This phenomenon has 
been observed in planetary atmospheres (Noll et 
al. 1988, Noll & Larson 1991, Fegley & Lodders 
1994, Bezard et al. 2002) and in brown dwarfs 
(Noll et al. 1997, Griffith & Yelle 1999, Saumon 
et al. 2000, 2006). In the convection zone that 
forms in the lower atmosphere, mixing is vigor- 
ous and occurs on a time scale T m ; x = H p /v conv 
where H p is the pressure scale height and v con v is 
the convective velocity, obtained from the mixing 
length theory. The physical process of mixing in 
the radiative zone is not known but it should take 
place on a much longer time scale. For lack of a 
theory of mixing in brown dwarf atmospheres, it is 
assumed to be a diffusive process and the mixing 
time scale is given by T m ; x = H p /K ZZl where K zz 
is the coefficient of eddy diffusivity, considered a 
free parameter here. For the present study, mod- 
els have been calculated with K zz = (no mixing, 
chemical equilibrium), 10 2 and 10 4 cm 2 s _1 in the 
radiative zone. For comparison, the mixing time 
scale in the convection zone would correspond to 
values of K zz w 10 7 -10 9 cm 2 s -1 . 
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We have generated a large grid of models with 
ranges of T e s, surface gravity g, f sc d, metallicity 
[M/H], and mixing coefficient K zz . In this paper 
we treat the dynamical parameters, K zz and / se di 
as independent variables, in the same way as we 
treat g and T e ff. However, the global properties of 
each dwarf must ultimately define the true prop- 
erties of the cloud, and the vigor of atmospheric 
mixing. Two otherwise identical L dwarfs, for ex- 
ample, could not exhibit two different sedimen- 
tation efficiencies and eddy diffusion coefficients. 
Eventually we will be able to treat K zz and / sc d 
as functions of other physical properties, but this 
cannot happen until we understand how differ- 
ences in metallicity, rotation rate, effective tem- 
perature, and gravity impact sedimentation and 
mixing in the atmosphere. Thus for the time be- 
ing we treat K zz and / so d as independent vari- 
ables, and use them to understand the ranges of 
behaviour manifested by the known objects. 

Saumon ct al. (2006) have validated the mod- 
els by reproducing the entire red through mid- 
infrared spectrum of the T7.5 dwarf Gl 570D. 
Cushing et al. (in preparation) also find very good 
agreement between the model spectra and the ob- 
served zJHKL 1 and IRS spectra of L0-T4 dwarfs. 
Work in progress comparing the calculated and 
observed near-infrared colors of L and T dwarfs 
shows that L dwarfs are well represented by mod- 
els with / sc d = 1-2, while T dwarf colors are well 
matched by models with / sc d = 4 or by cloudless 
models (i.e., the clouds lie below the photosphere 
at all wavelengths). Finally we note that because 
the models are complex and CPU-intensive, not 
all combinations of all parameters are available at 
this time. 

4.2. Comparison of Models and Data 

Figures 4, 5, and 6 each show four color-color 
diagrams that compare new and previously pub- 
lished 2.2-7.9 jiva photometry with synthetic col- 
ors generated from models covering appropriate 
ranges of T c ff, / sc d, log g, and K zz . Generally, the 
model sequences are generated for 500 < T c g < 
2400 K in steps of 100 K. However, the dusti- 
est (/ so d = 1 and 2) models with vertical mixing 
(K zz = 10 2 and 10 4 cm 2 s -1 ) are generated for the 
smaller range 800 < T e g < 2400 K because cal- 
culating cold, dusty, atmospheres with mixing is 
difficult and CPU intensive. This limitation does 



not affect our analysis or conclusions because the 
condensate cloud decks are expected to lie below 
the photosphere for T c g < 1000 K. Consequently, 
we do not consider inappropriate dusty models for 
mid- to late-T dwarfs or cloudless models for early- 
to mid-L dwarfs. All model sequences assume so- 
lar metallicity, except for one metal-poor cloudless 
sequence. 

Figure 4 demonstrates the colors' sensitivity to 
gravity for the two extremes of cloud properties: 
a model with / sc d = 1 and a completely cloud- 
less model. Both models assume no vertical mix- 
ing. The gravities shown (log g = 4.5, 5.0 and 
5.5) correspond to ~ 15-75 Jupiter masses for 
ages greater than 0.3 Gyr (Burrows et al. 1997), 
which are appropriate for our field sample with 
700 < T cS < 2200 K. Figure 5 shows sequences 
for log g = 5.0 and a full range of sedimention 
parameters, / sc d =1, 2, 3, 4, and cloudless, with- 
out vertical mixing. It also shows a metal-poor 
([M/H]= —0.3), logg = 5.0, cloudless sequence. 
Figure 6 shows / sc d =1, 2, 3, and cloudless se- 
quences for \ogg = 5.0, with vertical mixing pa- 
rameters of K zz = 10 2 and 10 4 cm 2 s _1 . 

In general, the varying colors of brown-dwarf 
model spectra are not straightforward functions of 
T e g, gravity, composition, and / sc d. One must also 
consider the variations in chemistry and overall 
gas and cloud opacity that affect the depth of the 
"photosphere." The relative strengths of absorp- 
tion bands in the spectrum are complex functions 
of all these parameters. For the sake of brevity 
and clarity, we will explain the trends for those 
cases where the interpretation is clear. 

4-2.1. Gravity 

Figure 4 indicates that gravity only impacts 
significantly (at the ^20% level) the K-L', K- 
[3.55], and [4.49]-[5.73] colors of mid- to late-T 
dwarfs. The [5.73]-[7.87] color is not significantly 
impacted, nor are any of the colors of dusty or 
cloudless dwarfs earlier than type T6 (although 
the absolute magnitudes can differ, as we discuss 
in §5). 
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Fig. 4. — Color-color diagrams and the effect of gravity: (a) K-L' vs. L'-M', (b) K-[3.55] vs. [3.55]-[4.49], (c) [3.55]-[4.49] vs. 
[4.49]— [5.73], and (d) [4.49]— [5.73] vs. [5.73]— [7.87]. Circles represent M dwarfs, triangles represent L dwarfs, and squares represent 
T dwarfs. These dwarfs are identified: 2MASS J0937347+293142, 2MASS J12171110-0311131, 2MASS J2244316+204343, SDSS 
J080531. 80+481233 and SDSS J151643.01+305344.4. The curves show the synthetic colors as a function of T cff for two cloud 
models (/ se< j = 1 - orange curves; cloudless - black curves), three gravities (logg = 4.5 - short-dashed curves; log 3 = 5.0 — 
solid curves; logg = 5.5 — long-dashed curves), and fixed solar metallicity. In panels (a)-(c), T e g = 2400 K lies at the bottom 
and T c ff decreases upward. In panel (d) the color sequence progresses from left to right as T a g decreases. 
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For T ff 1000 K, the K-b&nd flux is sensitive 
to gravity because it corresponds to a peak in the 
H2 CIA opacity. This opacity depends quadrati- 
cally on the gas number density, unlike most other 
opacity sources that are roughly linear functions of 
the density. Higher gravities cause higher photo- 
spheric pressures and a greater CIA opacity, which 
depresses the if -band flux. Figures 4a and 4b 
indeed show redder K-L' and iiT— [3.55] colors as 
log g increases. Only the coolest T dwarfs are af- 
fected, however, because at higher T e s the H 2 
opacity in the i-T-band increases and CIA is no 
longer dominant. 

The effect of gravity on the IRAC colors (Fig- 
ures 4c and 4d) is more subtle, and the cause varies 
with T e g or spectral type. The chemistry of CH4 
and PH 3 , and the temperature-pressure profile, 
are all dependent on gravity. The net result is 
bluer [4. 49] -[5. 73] colors as gravity increases, for 
T off <1400 K, or late-L and all T types. 

4-2.2. Sedimentation efficiency 

Figure 5 shows that the sedimentation effi- 
ciency, which regulates the thickness of the cloud 
decks and the geometric depth of the photosphere, 
has a large effect on the near- and mid-infrared 
colors of ultracool dwarfs. The wavelength de- 
pendence of the cloud opacity is fairly flat in the 
infrared (Ackerman & Marley 2001), except for 
a possible silicate feature around 10 /xm (Cushing 
et al. 2006). Hence, the cloud affects all bands 
from the optical to the mid-infrarcd. At longer 
wavelengths, the gas opacity tends to dominate 
and the mid-infrared flux is less sensitive to the 
choice of / se d, except when the cloud is very thick 

(/sod » 1). 

The range of / so d depicted in Figure 5 lies 
between the cloudless and thick-cloud extremes 
shown in Figure 4, and so the model sequences 
have intermediate colors. All the model color 
sequences in Figure 5 show a sharp reversal at 
T c ff « 1300 K. This reversal is also seen in the 
observed colors at the L-T transition (triangle 
to square symbols in the plots), where the atmo- 
spheres are changing rapidly from cloudy to cloud- 
less. The thick-cloud model (/ so d = 1) generally 
gives the closest fit to all the color-color diagrams, 
except for K-L' vs. L'-M' , which is better fit with 
cloudless models. Nonetheless, a significant mis- 
match remains, which compels us to explore model 



parameters other than gravity and /sed, in search 
of better agreement with the observations. 

4.2.3. Metallicity 

Figure 5 shows that decreasing the metallicity 
of the cloudless models by 0.3 dex impacts only 
the colors of the latest T dwarfs (T cff £ 1000 K) 
between 2.2 and 3.8 /im. Reducing the metallic- 
ity has a qualitatively similar effect on the chem- 
istry of C, N, and O as increasing the gravity or 
gas density (Lodders & Fegley 2002). The abun- 
dances of molecules composed of two metals (e.g. 
CO) decrease much faster than the abundances of 
molecules with only one metal (e.g. H2O). Thus, 
CO and N 2 will be relatively less abundant than 
H2O, CH4 and NH3. In addition, a decrease in 
the metallicity reduces the H2O abundance and 
the overall opacity in the atmosphere (which is 
dominated by H 2 absorption), and results in in- 
creased gas density. Finally, the importance of H2 
CIA increases in a metal poor atmosphere. All 
these effects mimic the effect of increasing gravity. 
The [M/H]= —0.3, logg = 5 sequence (Figure 5) is 
very similar to the [M/H]=0, logg = 5.5 sequence 
(Figure 4) in all the color-color diagrams shown. 
In a future paper, we will explore whether de- 
creasing metallicity or increasing gravity in dusty 
L dwarfs has effects similar to those of increasing 
sedimentation efficiency (i.e., thinning of the cloud 
decks), and vice versa. 

4.2.4. Vertical mixing 

Figures 4 and 5 show that the models with rea- 
sonable ranges of logg, / sc d, and [M/H] do not 
reproduce the observed colors of L and T dwarfs 
very well, although the general trends with T c g are 
satisfactory. The worst mismatches are seen in the 
K-L' and/or L'-M' colors (Figures 4a and 5a), 
where discrepancies of up to a full magnitude are 
seen for the T dwarfs. The synthetic [5.73]-[7.87] 
colors (Figures 4d and 5d) are also ~ 0.3 mag too 
blue. This discrepancy is consistent with the de- 
tailed fit of the SED of the T7.5 dwarf Gl 570D 
with the same models (Saumon et al. 2006) where 
the fluxes in the deep CH4 and H2O bands from 
5.5 to 8 /im (Figure 1) are underestimated. 
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Fig. 5. — Color-color diagrams and the effect of sedimentation efficiency. The diagrams and symbols are as 
in Figure 4. The curves show the synthetic colors as a function of T e g for five cloud models (/ se d = 4 - blue 
curves; / se d = 3 - green curves; / sc d = 2 - red curves; / sc d = 1 - orange curves; cloudless - black curves), 
two metallicities (solar - all solid curves; [M/H]— —0.3 - dashed curve), and fixed logg = 5.0. 
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It is also puzzling that the wide ranges of pa- 
rameters considered so far are unable to repro- 
duce the range in the observed colors of L dwarfs, 
shown as triangles in Figures 4 and 5. For exam- 
ple, the observed range in K —V is ~ 1 magnitude, 
and in [4.49] - [5.73] it is ~ 0.5 magnitudes (Fig- 
ure 3). The synthetic colors for the appropriate 
range of T cff (1400 K - 2400 K, G04) only vary by 
~ 0.6 magnitudes at K — V and by ~ 0.3 magni- 
tudes at [4.49] - [5.73], for 4.5 < log g < 5.5 and 

1 < /sed < 4. 

As described in §4.1, there is growing evidence 
that non-equilibrium chemistry caused by verti- 
cal mixing is common in T dwarfs. The spec- 
trum of the T6 dwarf Gl 229B shows enhanced 
CO and depleted NH 3 (Noll ct al. 1997; Oppen- 
heimer et al. 1998; Griffith & Yelle 1999; Saumon 
et al. 2000), and the mid-infrared spectrum of 
the T7.5 dwarf Gl 570D also shows depleted NH 3 
(Saumon et al. 2006). Models incorporating non- 
equilibrium chemistry have been shown to repro- 
duce the M' photometry of brown dwarfs bet- 
ter than equilibrium models (Leggett et al. 2002; 
G04). Another example is shown in Figure 1: 
the 3 /im CH4 band in the spectrum of the L9 
dwarf 2MASS J09083803+5032088 is much weaker 
than predicted by chemical equilibrium models 
with parameters appropriate for late-L field dwarfs 
(T off = 1400 K, / sod = 2, and log.? = 5). 

Figure 6 shows non-equilibrium models com- 
puted with two mixing efficiencies, K zz = 10 2 
and 10 4 cm 2 s -1 . The synthetic colors are very 
sensitive to vertical mixing, because the filter 
bandpasses sample absorption features caused by 
molecules whose abundances are significantly af- 
fected by the mixing. Figure 1 shows that K, V, 
and IRAC channels 1 and 4 sample CH4 bands; 
K, M', and IRAC channel 2 sample CO bands; 
and all filters (especially IRAC channels 2 and 3) 
are affected by H2O bands (Saumon et al. 2003). 
The red part of the IRAC channel 4 flux is also 
affected by NH3 absorption in the T dwarfs. Non- 
equilibrium effects on the CO and CH4 bands 
appear below T c s ~ 1600 K (or spectral types 
later than L5), at which point carbon is no longer 
fully sequestrated in CO. All four diagrams in 
Figure 6 show that the 2-8 [im photometry of 
late-L and T dwarfs is well reproduced by mod- 
els that include vertical mixing, especially when 
T e g ^ 1400 K (or spectral types later than about 



L7) . Most of the data are matched by models with 
K zz ~ 10 4 cm 2 s" 1 , with residuals of < 20%. Fur- 
thermore, the observed range in color for L and T 
types can be accommodated by variations in K zz . 

The best match to the K-V vs. L'—M' diagram 
(Figure 6a) is obtained with the model parame- 
ters K zz = 10 4 cm 2 s" 1 , logy = 5, [M/H]=0, and 
/scd > 2 (for T cff > 1400 K) or no clouds (for 
T c ff < 1400 K). The same parameters apply to 
the #-[3.55] vs. [3.55]-[4.49] diagram (Figure 6b), 
except that / sc d = 2 provides the best fit over 
the entire range of T c g. Higher values of / se d give 
equally good agreement if the gravity is increased. 
The models are less successful in reproducing the 
[3.55]-[4.49] vs. [4.49]-[5.73] diagram (Figure 6c). 
Vertical mixing has little effect in this diagram, 
and /scd ~ 1-2 best matches the colors of the M 
and L dwarfs. Further improvement at the warm 
end of the color sequence can be achieved with 
above-solar metallicity at the expense of a worse fit 
at T-dwarf temperatures. Finally, the [4.49]-[5.73] 
vs. [5.73]-[7.87] diagram (Figure 6d) shows that 
the observed colors of the hotter objects ([4.49]- 
[5.73] <; 0.1) are best reproduced with / sc d = 1 
and K zz = 10 4 cm 2 s" 1 . 

To summarize, good agreement with all four 
color-color diagrams is obtained with models with 
solar metallicity, logg « 5, / sc d =1-2 for T e g > 
1400 K, and cloudless models at lower T e g, and an 
eddy diffusion coefficient of K zz = 10 4 cm 2 s _1 . 
The data strongly support the idea that non- 
equilibrium chemistry, as modelled by vertical 
transport, is prevalent in late-L and T dwarfs. 

4-2.5. Objects with unusual colors 

Several objects stand out in Figure 3 by being 
unusually blue or red for their spectral type in 
one or more colors. We use the model sequences 
in Figures 4-6 to make a preliminary assessment 
of these peculiar colors. 

• The L7.5 dwarf 2MASS J22443167+2043433 
is very red in K-V, V-M', #-[3.55], 
and [4.49]-[5.73], suggesting that its at- 
mosphere is very cloudy (/ se d < 1) and 
strongly affected by vertical mixing (K zz > 
10 4 cm 2 s" 1 ), as shown for example by Fig- 
ure 6(b). 
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Fig. 6. — Color-color diagrams and the effect of vertical mixing. The diagrams and symbols are as in 
Figure 4. The curves show the synthetic colors as a function of T c g for four cloud models (/ se d = 3 - 
green curves; / scc i = 2 - red curves; / sc d = 1 - orange curves; cloudless - black curves), two eddy diffusion 
coefficients {K zz = 10 2 cm 2 s _1 - solid curves; K zz = 10 4 cm 2 s _1 - dashed curves), solar metallicity, and 
logff = 5.0. 
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• SDSS J133148. 90-011651.4 is an L8 dwarf 
that is blue in [3.55]-[4.49] and [4.49]-[5.73], 
suggesting that it has a high value of K zz 
and thin clouds with / sc( j ^ 3 (Figure 6(b)). 

• SDSS J080531.80+481233 (L9.5) is blue in 
[4.49]-[5.73] and K-L', and is similar to 
SDSS J133148. 90-011651.4. It is best repro- 
duced by a high sedimentation model, with 
/sed ~ 4 or perhaps no clouds, and modest 
vertical mixing. 

• The T0.5 dwarf SDSS J151643. 01+305344.4 
is red in #"-[3.55] and [5.73]-[7.87], which 
requires both vertical mixing and fairly thick 
clouds with /sod < 2, as shown for example 
by Figure 6(d). 

• The well-known peculiar T6 dwarf 2MASS 
J09373487+2931409 has such extremely red 
K-L' and if-[3.55] colors that it can only be 
explained with a high gravity and low metal- 
licity atmosphere, plus vertical mixing (e.g. 
Figures 4(b), 5(b) and 6(b)). 

• 2MASS J12171110-0311131 (T7.5) is very 
blue in [3. 5 5] -[4. 49], moderately red in 
[4.49]-[5.73], and moderately blue in other 
colors (P06). This dwarf has a bright K- 
band flux peak, implying weak H 2 CIA 
and therefore low gravity (Burgasser et al. 
2006a). However, the color sequences in 
Figures 4 and 5 show that decreasing grav- 
ity and increasing metallicity would not re- 
produce the observed colors, so the high K 
flux may be due to something other than 
low H2 opacity. The colors may be due 
to extremely efficient vertical mixing, e.g., 
K zz ~ 10 6 cm 2 s" 1 , which leads to reduced 
CH4 and brighter K and [3.55] magnitudes 
(see Figure 6b). This T dwarf has been 
imaged by HST and appears to be single 
(Burgasser et al. 2006c), so duplicity does 
not seem to explain its unusual colors. 

As we note in §4.1, the above discussion im- 
plies that all five model parameters, T e g, grav- 
ity, [M/H], / sc d and K zz , can be individually var- 
ied to fit the properties of each unusual object. 
In fact the dynamical parameters / se d and K zz 
will depend on the other three bulk properties, 



as well as on rotation and the detailed chemi- 
cal composition. Further modeling of individual 
objects, such as 2MASS J22443167+2043433 and 
SDSS J133148.90-011651.4 which seemingly have 
similar spectral types and vertical mixing but very 
different / sc d, will help to illuminate such issues. 

4.3. Effects of Gravity and Vertical Mixing 
on Luminosity 

P06 suggest that the spread in the IRAC pho- 
tometry of T dwarfs is due to variations in gravity 
within the sample. Except when they are very 
young, brown dwarfs cool at nearly constant radii 
and their gravity strongly correlates with their 
mass and their luminosity. The effect of gravity 
on models is therefore better revealed in color- 
magnitude diagrams than in color-color diagrams. 

Figure 7 shows the IRAC channel 1 and 2 ab- 
solute magnitudes as functions of color, Mr 3 . 55 i vs. 
[3.55]-[4.49] and M [4A9] vs. [4.49]-[5.73], for those 
objects which have measured parallaxes. The data 
reveal clear sequences as a function of spectral 
type, with increased scatter in the T dwarf regime. 
Evolutionary calculations indicate that for a pop- 
ulation of very-low-mass stars and brown dwarfs 
older than iJ 1 Gyr (such as our present sample) , 
objects with lower T c g will span a wider range of 
mass and gravity (see Figure 1 of Saumon et al. 
2006). However, the model sequences shown in 
Figure 7 indicate that, for the plausible range of 
log g = 4.5 — 5.5, the width due to gravity of the 
modelled sequence is narrower than the observed 
T dwarf sequence, at least for the (M[4. 49 ]) vs. 
[4.49]-[5.73] diagram. We conclude that a range 
of gravities is not sufficient by itself to explain 
the width of the T dwarf sequence in these dia- 
grams. Increasing the vertical mixing coefficient 
from K zz = to 10 4 cm 2 s _1 shifts the T dwarfs' 
[3.55]-[4.49] colors blueward and the [4.5]-[5.8] 
colors redward, by several tenths of a magnitude 
- and mostly recovers the observed sequence. 
This large shift clearly shows that a range of ver- 
tical mixing efficiency, (K zz ~ 10 2 -10 6 cm 2 s _1 ) 
in cloudless brown-dwarf atmospheres could repro- 
duce the observed range of T dwarf 3-8 /zm colors. 
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Fig. 7. — IRAC color-magnitude diagrams for those objects with measured parallax. The symbols are 
described in Figure 4. Known binaries (with combined photometry) are labelled, as is the peculiar T6 dwarf 
2MASS J0937347+293142. Cloudy models with / sod = 2 and 1400 < T cS < 2400 K are shown by the 
red curves, and cloudless models with 500 < T e g < 1400 K are shown by the black curves. The plausible 
range of gravity for field brown dwarfs (log 3 = 4.5-5.5) is bracketed by the short-dashed and long-dashed 
curves, respectively. These sequences do not include vertical mixing. Models including vertical mixing, with 
K zz = 10 4 cm 2 s^ 1 and logg = 5, are shown by the solid curves. Labels on the right axes indicate T e g (K) 
for the non-equilibrium models. 
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5. Effective Temperature Indicators 



Of the colors considered in this paper, two show 
promise as effective temperature indicators. Fig- 
ure 3 shows that both K-L' and [3.55]-[4.49] ex- 
hibit relatively little scatter and vary monotoni- 
cally with spectral type. The monotonicity of K— 
L' with spectral type was first noted by Stephens 
et al. (2001). This relation was subsequently used 
by Chauvin et al. (2004) and Neuhauser et al. 
(2005) to assign a spectral type and T c g to candi- 
date planetary-mass companions of more massive 
objects. 

Figure 8 shows the observed K-L' and [3.55]- 
[4.49] colors as a function of T e ff, as determined 
by G04 and Luhman et al. (2006). Cloudy 
(/sod =2) and cloudless model color sequences are 
also shown. The observed colors are well repro- 
duced by the models in the appropriate range of 
T e ff when vertical mixing is included. Although 
[3.55]-[4.49] is constant down to ~ 1700 K, the 
onset of CH4 absorption moves the colors steadily 
redward through the late-L and the T sequences. 
As shown in the figure, the putative Y dwarfs (i.e., 
dwarfs cooler than T dwarfs) will also likely be red 
in both colors. 

We note that care must be taken when applying 
empirical color-T e ff relations appropriate for field 
objects to young low-mass and low-gravity, and 
possibly non-solar metallicity, objects. For late-T 
dwarfs, H2 CIA opacity is important, particularly 
in the K band, and extremely sensitive to gravity 
and metallicity (§4.2.1, 4.2.3). For example, for a 
mid-T dwarf with K-L' w 2.5, varying log g by 0.5 
dex produces a change in T e ff of 100 K, and vary- 
ing [M/H] by 0.3 dex a change of ^50 K. In ad- 
dition, K-L' is sensitive to the mixing parameter 
K zz for late-L to early-T types, and [3.55]-[4.49] is 
also sensitive to this parameter for late-L to late-T 
types. Figure 8 shows that if K zz is substantially 
larger or smaller in the atmospheres of low-mass 
objects, errors in derived T c g of 100 — 300 K will 
result. 

6. Conclusions 

Our NIRI V (3.75 /mi), W (4.70 /tm), and 
IRAC channels 1-4 (3.55, 4.49, 5.73 and 7.87 /xm) 
photometry, when combined with previously pub- 
lished data, constitutes a sample of 30-50 objects 
(depending on the colors) that is large enough to 




3000 2500 2000 1500 1000 500 

T K 

eff 



Fig. 8. — Colors as indicators of T e ff. The values 
of T e ff are determined from bolometric luminosi- 
ties by G04 and Luhman et al. (2006); the G04 
values have been revised for the binarity of Kelu-1 
(Liu & Leggett 2005), LHS 102BC (Golimowski et 
al. 2004b), SDSS J042348. 57-041403.5 (Burgasser 
et al. 2005) and SDSS J102109.69-030420.1 (Bur- 
gasser et al. 2006c). The symbols are described 
in Figure 4. Approximate spectral types are indi- 
cated, using the T c ff vs. spectral type relationship 
of G04. The errors in T e g are dominated by un- 
certainty in age. The curves show model color 
sequences with [M/H] = 0, log<? = 5, and cloud- 
less atmospheres with 500 < T G ff < 1400 K (black 
curves) or cloudy atmospheres with / sc d = 2 and 
1400 < T cff < 2400 K (red curves). The dashed 
curves are chemical equilibrium models (K zz = 0) 
and the solid curves are non-equilibrium models 
with K zz — 10 4 cm 2 s _1 . The peculiar T6 dwarf 
2MASS J0937347+293142 is identified. 
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establish the behaviour of the entire L-T spec- 
tral sequence in color-color diagrams. Useful com- 
parisons with synthetic mid-infrared colors are 
now possible. Current model atmospheres are 
quite complex, incorporating detailed molecular 
and atomic opacities, the formation and sedimen- 
tation of condensates, and vertical mixing leading 
to non-equilibrium chemistry. 

Generally the models reproduce the observed 
mid-infrared colors well. For early- to mid-L 
dwarfs (T cff > 1600 K), the synthetic IRAC colors 
are insensitive to variations in the model parame- 
ters, except for the case of thick clouds (/ se d = !)• 
The colors of the late-L and T dwarfs are only re- 
produced by models that include non-equilibrium 
chemistry arising from vertical transport (mixing) 
in the radiative region of the atmosphere. Effi- 
cient mixing is necessary, with an eddy diffusion 
coefficient of K zz « 10 4 cm 2 s _1 . The comparisons 
presented here show that vertical mixing in the at- 
mospheres of late-L and T dwarfs (T e s ^ 1600 K) 
affects the chemistry of carbon and is widespread 
among these objects. 

Late-L and early-T dwarfs that are unusually 
red or unusually blue in the near-infrared are also 
red or blue, respectively, in some of the mid- 
infrared colors. The red objects are best repro- 
duced by non-equilibrium models with very low 
sedimentation efficiency (/ so d ~ 1) and hence thick 
condensate cloud decks, and the blue objects are 
reproduced by non-equilibrium models with very 
high sedimentation efficiency (/ se d ~ 3~4) and 
hence thin condensate cloud decks. 

The models account for most of the observed 
range in color and absolute magnitude via vari- 
ations in the eddy mixing coefficient within the 
range K zz = 10 2 -10 6 cm 2 s" 1 . The mixing mech- 
anism in the radiative region of brown dwarf atmo- 
spheres is presently unknown, so it is not possible 
to assess whether such a wide range is reasonable 
or not. The remaining < 20% discrepancies be- 
tween observed and modelled colors and magni- 
tudes can be attributed to several possible causes, 
including variations in cloud properties that are 
not captured in our one-dimensional cloud model, 
uncertainties in the opacities of CH4 and NH 3 , the 
choices of C, N, and O abundances, and perhaps 
variations in the C/O ratio in brown dwarf atmo- 
spheres. 

Vertical mixing may also be an important con- 



sideration for the direct detection of giant plan- 
ets around nearby stars (G04; Marley et al. 2006, 
Hinz ct al. 2006). After the discovery of Gl 229B, 
Marley et al. (1996) suggested that a substantial 
4-5 jj,m flux peak should be a universal feature of 
giant planets and brown dwarfs. This expectation, 
combined with a favorable planet/star flux ratio, 
has made the M band a favorite for planet detec- 
tion (Burrows et al. 2005, Marley et al. 2006). 
However, Figures 4 and 8 show that the M or 
4.6 ixm region is fainter - and the L or 3.6 fim 
region is brighter - than the equilibrium chem- 
istry models predict, at least for T c g > 700 K. 
Given these and other considerations (such as 
cleaner atmospheric transmission at L compared 
with M) the comparative advantage of ground- 
based searches for young, bright giant planets at 
M may be somewhat less than currently expected 
(see also Marley et al. 2006). Another implica- 
tion of the reduced M-band flux is that planetary 
masses (and upper limits) inferred from M-band 
measurements using equilibrium chemistry models 
(Hinz et al. 2006) will be too low. 

This initial study is a comparison of observed 
and synthetic broadband color-color diagrams of 
brown dwarfs, and it is limited to identifying and 
interpreting trends in the data. The model pa- 
rameters that best match the colors of a specific 
object vary with color, which, unless the models 
are perfect, limits the level of analysis that can 
be performed. Being differential in nature, col- 
ors can magnify small errors in models that oth- 
erwise agree well with the data. For example, the 
<~ 0.3 mag mismatch between the observed and 
synthetic [5.73]-[7.78] colors of latc-T dwarfs (Fig- 
ures 4(d), 5(d), and 6(d)) betrays the generally 
good fit of the same models with the entire SED of 
the T7.5 dwarf Gl 570D (Saumon et al. 2006). Be- 
cause of the haphazard overlap of the photometric 
and molecular bands and the complex variations 
in brown dwarf spectra, color-color diagrams can 
be difficult to interpret. With global trends identi- 
fied and generally good agreement between mod- 
els and data, the detailed study of both typical 
and anomalous brown dwarfs can be tackled. Our 
analysis of ground-based and Spitzer photometry 
and spectra of such objects is underway. 

We are grateful to Malcolm Currie at the 
Rutherford Applcton Laboratory, and the Starlink 



17 



team, for help with the NIRI data analysis. We 
are also grateful to Linhua Jiang of the University 
of Arizona for help with the IRAC data analysis. 
Brian Patten of CfA Harvard was very generous 
with pre-publication access to IRAC data and in 
helping with the IRAC data analysis. 

This work is based in part on observations made 
with the Spitzer Space Telescope, which is oper- 
ated by the Jet Propulsion Laboratory Califor- 
nia Institute of Technology under a contract with 
NASA, and on observations made at the Gem- 
ini Observatory, which is operated by the Asso- 
ciation of Universities for Research in Astronomy 
Inc. (AURA), under a cooperative agreement with 
the NSF on behalf of the Gemini partnership: the 
National Science Foundation (United States), the 
Particle Physics and Astronomy Research Council 
(United Kingdom), the National Research Coun- 
cil (Canada), CONIC YT (Chile), the Australian 
Research Council (Australia), CNPq (Brazil) and 
CONICET (Argentina). 

Support for this work was provided by NASA 
through an award issued by JPL/Caltech. This 
work was also supported in part under the aus- 
pices of the U.S. Department of Energy at Los 
Alamos National Laboratory under Contract W- 
7405-ENG-36. MSM acknowledges the support of 
the NASA Office of Space Sciences. TRG is sup- 
ported by the Gemini Observatory, which is oper- 
ated by AURA on behalf of the international Gem- 
ini partnership listed above. SKL acknowledges 
the support of the Gemini Observatory, the Joint 
Astronomy Centre and the U.K. Particle Physics 
and Astronomy Research Council. 

Funding for the SDSS and SDSS-II has been 
provided by the Alfred P. Sloan Foundation, 
the Participating Institutions, the National Sci- 
ence Foundation, the U.S. Department of En- 
ergy, NASA, the Japanese Monbukagakusho, the 
Max Planck Society, and the Higher Education 
Funding Council for England. The SDSS Web 
Site is http://www.sdss.org/. The SDSS is man- 
aged by the Astrophysical Research Consortium 
for the participating institutions: American Mu- 
seum of Natural History, Astrophysical Institute 
Potsdam, University of Basel, Cambridge Uni- 
versity, Case Western Reserve University, Univer- 
sity of Chicago, Drexel University, Fermilab, In- 
stitute for Advanced Study, Japan Participation 
Group, Johns Hopkins University, Joint Institute 



for Nuclear Astrophysics, Kavli Institute for Parti- 
cle Astrophysics and Cosmology, Korean Scientist 
Group, Chinese Academy of Sciences (LAMOST), 
Los Alamos National Laboratory, Max-Planck- 
Institutc for Astronomy (MPIA), Max-Planck- 
Institute for Astrophysics (MPA), New Mexico 
State University, Ohio State University, University 
of Pittsburgh, University of Portsmouth, Prince- 
ton University, United States Naval Observatory, 
and University of Washington. 

References 

Ackerman, A. S. & Marley, M. S. 2001, ApJ, 556, 
872 

Beichman, C. A., Chester, T. J., Skrutskie, M., 
Low, F. J. & Gillett, F. 1998, PASP, 110, 480 

Bezard, B., Lellouch, E., Strobcl, D., Maillard J.- 
P. & Drossard, P. 2002, Icarus, 159, 95 

Burgasser, A. J. et al. 2002, ApJ, 564, 421 

Burgasser, A. J, Marley, M. S., Ackerman, A. S., 
Saumon, D., Lodders, K., Dahn, C. C, Harris, 
H. C, Kirkpatrick, J. D. 2002, ApJ571, L151 

Burgasser, A. J., Reid, I. N., Leggett, S. K., Kirk- 
patrick, J. D., Liebert, J., Burrows, A. 2005, 
ApJ, 634, 177 

Burgasser, A. J., Burrows, A., Kirkpatrick, J. D. 
2006a, ApJ, 639, 1095 

Burgasser, A. J., Geballe, T. R., Leggett, S. 
K., Kirkpatrick, J. Davy & Golimowski, D. A. 
2006b, ApJ, 637, 1067 

Burgasser, A. J., Kirkpatrick, J. D., Cruz, K. L., 
Reid, I. N., Leggett, S. K., Liebert, J., Burrows, 
A., Brown, M. E. 2006c, ApJ, in press, astro- 
ph/0605577 

Burrows, A., et al. 1997, ApJ, 491, 856 

Burrows, A., Hubbard, W. B., Lunine, J. I., 
Liebert, J. 2001, Reviews of Modern Physics, 
73, 719 

Burrows, A., Sudarsky, D., Hubcny, I. 2006, ApJ, 
640, 1063 

Cavanagh, B., Hirst, P., Jcnncss, T., Economou, 
F., Currie, M. J., Todd, S. & Ryder, S. D., 2003, 
Astronomical Data Analysis Software and Sys- 
tems XII, ASP Conference Series, Vol. 295, H. 
E. Payne, R. I. Jedrzcjcwski, and R. N. Hook, 



18 



eds., 237 

Chauvin, G., Lagrange, A.-M., Dumas, C, Zuck- 
erman, B., Mouillet, D., Song, I., Beuzit, J.-L., 
& Lowrance, P. 2004, A&A, 425, L29 

Chiu, K., Fan, X., Leggett, S. K., Golimowski, D. 
A., Zheng, W., Geballe, T. R., Schneider, D. P. 
& Brinkmann, J. 2006, AJ, 131, 2722, (C06) 

Cushing, M. C. ct al. 2006, ApJ, 648,614 

Dahn C C. et al. 2002, AJ, 124, 1170 

Fazio, G. et al., 2004, ApJS, 154, 10 

Fegley, B. Jr. & Lodders, K. 1994, Icarus, 110, 
117 

Fegley, B. Jr. & Lodders, K. 1996, ApJ, 472, L37 

Geballe, T. R. et al. 2002, ApJ, 564, 466 

Golimowski D. A. et al. 2004a, AJ, 127, 3516, 
(G04) 

Golimowski D. A. et al. 2004b, AJ, 128, 1733 

Griffith, C.A. & Ycllc, R.V. 1999, ApJ, 519, L85 

Hinz, P.H., Heinze, A.N., Sivanandam, S., 
Miller, D.L., Kenworthy, M.A., Brusa, G., Freed, 
M. & Angel, J.R.P. 2006, ApJ, in press, astro- 
ph/0606129 

Hodapp, K. W. et al. 2003, PASP, 115, 1388 

Houck, J. et al. 2004, ApJS, 154, 18 

Kendall, T. R., Delfossc, X., Martin, E. L., 
Forvcillc, T. 2004, A&A, 416, L17 

Knapp, G. R. ct al. 2004, AJ, 127, 3553, (K04) 

Leggett, S. K. et al. 2002, ApJ, 564, 452 

Leggett, S. K., Hawarden, T. G., Currie, M. J., 
Adamson, A. J., Carroll, T. C, Kerr, T. H., 
Kuhn, O. P., Seigar, M. S., Varricatt, W. P. & 
Wold, T. 2003, MNRAS, 345, 144 

Liu, M. C. & Leggett, S. K. 2005, ApJ, 634, 616 

Lodders, K. & Fegley, B., Jr. 2002, Icarus, 393, 
424 

Luhman, K. L. et al. 2006, ApJ, in press 

Marley, M.S., Fortney, J., Seager, S., Barman, T. 
2006, in Protostars and Planets V, B. Reipurth, 
Ed., in press. 

Marley, M. S., Saumon, D., Guillot, T., Freedman, 
R. S., Hubbard, W. B., Burrows, A., & Luninc, 



J. I. 1996, Science, 272, 1919 

Marley, M. S., Seager, S., Saumon, D., Lodders, 
K., Ackerman, A. S., Freedman, R. S. & Fan, 
X. 2002, ApJ, 568, 335 

Neuhauser, R., Guenther, E. W., Wuchterl, G., 
Mugrauer, M., Bcdalov, A., & Hauschildt, P. H. 
2005, A&A, 435, L13 

Noll, K. S., Geballe, T. R., Leggett, S. K. & Mar- 
ley, M. S. 2000, ApJ, 541, L75 

Noll, K.S., Geballe, T.R. & Marley, M.S. 1997, 
ApJ, 489, 87 

Noll, K.S., Knacke, R.F., Geballe, T.R. & Toku- 
naga, A.T. 1988, ApJ, 324, 1210 

Noll, K.S., & Larson, H.P. 1991, Icarus, 89, 168 

Patten, B. M., et al. 2006, ApJ, in press, astro- 
ph/0606432 

Rcid, I. N. & Cruz, K. L. 2002, AJ, 123, 466 

Rocllig, T. L., ct al. 2004, ApJS, 154, 418 

Saumon, D., Geballe, T. R., Leggett, S. K., Mar- 
ley, M.S., Freedman, R. S., Lodders, K., Fegley, 
B., Jr. & Sengupta, S. K. 2000, ApJ, 541, 374 

Saumon, D., Marley, M.S., Cushing, M.C., Leggett, 
S.K., Rocllig, T.L., Lodders, K. & Freedman, 
R.S. 2006, ApJ, in press, astro-ph/0605563 

Saumon, D., Marley, M. S., Lodders, K. & Freed- 
man, R. S. 2003, Brown Dwarfs, Proceedings of 
IAU Symposium 211, ed. E. Martin. San Fran- 
cisco: Astronomical Society of the Pacific, 2003, 
345 

Simons, D. A. & Tokunaga, A. 2002, PASP, 114, 
169 

Skrutskie, M. F. et al. 2006, AJ, 131, 1163 

Stephens, D. C, Marley, M. S., Noll, K. S., & 
Chanovcr, N. 2001, ApJ, 556, L97 

Tokunaga, A. T., Simons, D. A. & Vacca, W. D. 
2002, PASP, 114, 180 

Toon, O.B., McKay, CP., Ackerman, T.P. & San- 
thanam, K. 1989, JGR 94, 16287 

Vrba, F. J. et al. 2004, AJ, 127, 2948 

Werner, M., et al. 2004, ApJS, 154, 1 

York, D. G. et al. 2000, AJ, 120, 1579 



19 



Table 1 

Gemini NIRI V and M' Photometry (MKO System) 



Name 


Type 1 


L (error) 


Date 


M (error) 


Date 


Notes 






(mag) 


(YYMMDD) 


(mag) 


(YYMMDD) 




SDSS J000013.54+255418.6 


T4.5 


13.03(0.03) 


040826 


13.28(0.10) 


040801 




2MASS J00345157+0523050 


T6.5 


13.30(0.04) 


040827 


12.66(0.10) 


040827 




2MASS J02431371-2453298 


T6 






13.27(0.16) 


050128 




SDSS J080531. 80+481233 


L9.5 


12.34(0.04) 


041230 


12.81(0.07) 


041230,050115 


2,3 


SDSS J083008. 12+482847.4 


L9 




12.75(0.15) 


040405 




SDSS J085758. 45+570851. 4 


L8 






11.55(0.10) 


040405 


4,5 


SDSS J111010.01+011613.1 


T5.5 


13.89(0.05) 


050115 






4 


SDSS J133148.90-011651.4 


L8 


12.73(0.03) 


050115 






2 


2MASS J15031961+2525196 


T5 






11.92(0.04) 


040509 


6 


2MASS J15530228+1532369AB 


T7 


13.83(0.05) 


050115 


13.29(0.15) 


040406 


7 


2MASS J16322911 + 1904407 


L7.5 






13.07(0.15) 


040530 




2MASS J22443167+2043433 


L7.5 






11.90(0.03) 


040801 


4 



1 Spectral types are based on near- infrared classification schemes for L dwarfs (Geballe et al. 2002) and T dwarfs (Burgasscr 
ct al. 2006b). 

2 Near-infrared colors are unusually blue for this spectral type. 

3 G04 reported V = 12.31 ± 0.05. 

4 Near-infrared colors are unusually red for this spectral type. 

5 Leggett et al. (2002) reported M' = 11.50 ± 0.05. 

6 G04 reported M 1 = 12.25 ± 0.15. 

7 T7 type applies to unresolved pair. Burgasscr et al. (2006c) report T6.5+T7. 



Table 2 

Filter Bandpasses at Operational Temperatures 



Filter 


Nominal 


50% Cut- 


50% Cut- 


Name 


Wavelength (pm) 


On (fim) 


Off (/mi) 


^MKO 


3.75 


3.43 


4.11 


"mko 


4.70 


4.57 


4.80 


IRAC channel 1 


3.55 


3.18 


3.92 


IRAC channel 2 


4.49 


4.00 


5.01 


IRAC channel 3 


5.73 


5.02 


6.43 


IRAC channel 4 


7.87 


6.45 


9.33 
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Table 3 
Spitzer IRAC Photometry 



Name 


Type 1 




Channel 1 






Channel 2 






Channel 3 






Channel 4 




Date 
(YYMMDD) 


Notes 


mjy 


[3.55] 


% err 


mjy 


[4.49] 


% err 


mjy 


[5.73] 


7o err 


mjy 


[7.87] % 


err 


SDSS J000013. 54+255418. 6 


T4.5 


0.91 


13.72 


1.0 


1.06 


13.07 


2.0 


1.09 


12.56 


10.0 


0.64 


12.50 


1.0 


050725 


2 


SDSS J075840. 33+324723.4 


T2 


3.05 


12.41 


0.4 


2.52 


12.13 


0.5 


2.37 


11.71 


1.8 


1.89 


11.33 


4.0 


051023 




SDSS J080531. 80+481233 


L9.5 


2.98 


12.44 


0.5 


1.92 


12.43 


0.5 


1.36 


12.32 


1.1 


0.93 


12.10 


1.8 


051022 


3 


SDSS J105213. 51+442255. 7 


T0.5 


1.43 


13.23 


0.1 


1.03 


13.10 


0.2 


0.98 


12.67 


0.5 


0.73 


12.36 


0.8 


051127 




SDSS J115553. 86+055957. 5 


L7.5 


2.07 


12.83 


0.2 


1.33 


12.83 


0.3 


1.15 


12.50 


1.0 


0.80 


12.26 


2.0 


051229 




SDSS J120747. 17+024424.8 


TO 


2.01 


12.86 


0.2 


1.46 


12.73 


0.3 


1.24 


12.42 


1.5 


0.82 


12.23 


5.0 


060706 




SDSS J151643. 01+305344. 4 


TO. 5 


1.01 


13.61 


0.5 


0.79 


13.39 


1.0 


0.80 


12.89 


1.0 


0.66 


12.47 


2.0 


050724 


4 


SDSS J152039. 82+354619. 8 


TO 


1.87 


12.94 


0.5 


1.25 


12.89 


0.2 


1.16 


12.49 


1.2 


0.80 


12.26 


0.8 


050726 




2MASS J22443167+2043433 


L7.5 


3.22 


12.35 


0.4 


2.58 


12.11 


0.5 


2.65 


11.59 


0.9 


1.83 


11.36 


1.5 


050723 


5 



1 Spectral types are based on near- infrared classification schemes for L dwarfs (Geballe et al. 2002) and T dwarfs (Burgasser et al. 2006b). 

2 Nearby bright star caused more uncertain photometry. Sky values were interpolated between the diffraction spikes of the bright star. 

3 Near-infrared colors are unusually blue for this spectral type. 

4 Near-infrared colors are unusually red for this spectral type. 

5 Near-infrared colors are extremely red for this spectral type. 



